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We build on the study of single-photon heat conduction in electronic circuits taking into account 
the back-action of the superconductor-insulator-normal-metal thermometers. In addition, we show 
that placing capacitors, resistors, and superconducting quantum interference devices (SQUIDs) into 
a microwave cavity can severely distort the spatial current profile which, in general, should be 
accounted for in circuit design. The introduction of SQUIDs also allows for in situ tuning of the 
photonic power transfer which could be utilized in experiments on superconducting quantum bits. 



I. INTRODUCTION 

In the past decade great progress has been achieved 
in the fields of circuit quantum electrodynamics (1-8) 
and photonic heat conduction (9-13). The first attempts 
to combine them by studying single-photon heat con- 
duction in a microwave cavity were recently reported in 
Ref. 14. This approach has the potential to deliver prac- 
tical benefits, both for investigating fundamental quan- 
tum phenomena and in order to deliver the practical tools 
needed in the field of superconducting quantum comput- 
ing where the manipulation of microwave photons is be- 
coming increasingly important 

In Ref. 14, a model to study single photons in a well- 
defined environment was introduced. The practical setup 
proposed allows for heat exchange between two normal- 
metal resistors in a superconducting cavity, for which the 
photonic heat conduction dominates. This scheme in- 
volves a pair of superconductor-insulator-normal-metal 
(SIN) tunnel junctions^^ to control the temperature of 
one of the resistors, enabling remote heating or cooling 
of low-temperature circuit elements. The coupling of the 
resistors to the superconducting cavity also allows them 
to act as engineered and controllable environments for 
the cavity modes. In practice, controlling and measuring 
the temperatures of the resistors inevitably introduces 
additional power sources for the electron clouds in the 
resistors. Additionally, there may also be heat transfer 
between the resistors by quasiparticle excitations in the 
superconductor. In this paper, we extend the analysis of 
the system in Ref. 14 by accounting for the SIN junc- 
tions and quasiparticle power sources and show that the 
observation of single-photon heating and cooling is still 
experimentally feasible. In addition, this more compre- 
hensive model allows us to impose a cutoff-temperature 
below which cooling is not observable. 

In the second half of this paper, we return to the clas- 
sical transmission line framework^^'^"* in order to inves- 
tigate how placing capacitors, resistors, and supercon- 
ducting quantum interference devices (SQUIDs) into the 
transmission line affects the photon mode profile. We 
show that the distortion of the photon modes responsi- 
ble for heat transport can be dramatic, and should poten- 
tially be accounted for in typical circuit design. Further- 
more, this mode distortion suggests a resolution to the 



discrepancy between the power transfer calculated with 
the quantum and classical models in the strong coupling 
limit, which was observed in Ref. 14. Since many future 
cavity photon experiments^^ will utilize two or more cav- 
ities, e.g., for filtering or isolating modes, we also focus 
on the effect of placing a capacitive break into the central 
superconducting strip of the transmission line. 

We study the possibility of using SQUIDs for in situ 
tuning of the single-photon power in a cavity. We fur- 
ther develop the classical transmission line model, and 
demonstrate that such tuning is indeed possible and can 
be rather effective. In the experiment of Ref. 9, in situ 
tuning of photonic heat conduction was demonstrated ex- 
perimentally by using SQUIDs to vary the impedance of 
the connecting circuit. However, no cavity was present 
and hence a lumped element model for the impedance 
was sufficient. In superconducting quantum computing, 
SQUIDs arc often placed into a cavity in order to cre- 
ate or readout quantum bits, qubits^''"'^'^ . Introducing 
SQUIDs into the line also provides some control over 
the eigenfrequencies of the photons in the cavity so that 
they may be, e.g, tuned in and out of resonance with a 
qubit. This tuning has been demonstrated experimen- 
tally through two approaches: either the SQUIDs termi- 
nate the transmission line, thus giving variable boundary 
conditions^®^'^^ , or they are placed in the center conduc- 
tor of the line''". In Ref. 28, the dynamical Casimir effect 
was demonstrated; photons were created in a SQUID- 
terminated transmission line by rapid variation of the 
SQUID flux. Tunable resonators are also required in sev- 
eral other practical applications: in Ref. 29, two photon 
number states are held in a cavity at different frequencies 
allowing for a parametric interaction between them. 

This paper is outlined as follows. In Sec. II, we review 
the basic theory required to calculate the heat power be- 
tween two normal-metal islands in the center conductor 
of the line and the theory of the SIN junctions needed 
to measure their temperature. We also discuss the heat 
power due to quasiparticles between two such islands into 
the model. We then study the remote heating and cooling 
in the presence of SIN junctions in Sec. III. Section IV 
gives a brief discussion of the discrete transmission line 
model and the mode profiles are solved with a capacitor 
placed into the line. In Sec. V, the model is expanded to 
include resistors and the power transfer between them is 
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FIG. 1. Schematic figure of two resistors, _Ri and R2, a ca- 
pacitor witfi capacitance D and a SQUID witli eflfective in- 
ductance Lj in a coplanar waveguide cavity. Tfie SIN probes 
are also indicated. Both resistors have a pair of thermometer 
probes and Ri has also a pair of refrigeration probes. This is 
an extension of the setup to study heat transfer in a cavity 
proposed in Ref. 14. 



rule with this Hamiltonian yields a photonic heat power 
between two such resistors of 
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where the first sum is over the cavity modes and the sec- 
ond over the photon number for each mode. The proba- 
bility that the /cth mode contains n photons is pj^j and the 
transition rate for photon number in the fcth mode due 



to the ith resistor is T^'^^^j^. For a transmission line of 
length L, inductance per unit length £, and capacitance 
per unit length c, the frequency of the fcth mode in the 
cavity is fk = k/{2L^/Ic), and the transition rates are 
given by 
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where 7^ = [2R,fim^{TTx,/ L)]/ {U), Xi is the position 
of resistor i = 1,2 and Ri its resistance. The function 
fB{u},T) = l/[exp(^) - 1] is the Bose-Einstein 
distribution function. 



calculated. In Sec. VI, we analyze the case of SQUIDs in 
the line by demonstrating the in-situ tuning of the cavity 
eigenfrequencies and power transfer. 



II. COMPARISON OF THE THERMAL POWER 
SOURCES 

A. Phononic and cavity photonic power 

In Ref. 14, the heat exchange between two normal- 
metal resistors in a superconducting cavity was studied 
in a setup similar to Fig. 1 with no SQUID or capac- 
itor in the line, i.e., Lj = and D — >■ 00. The equi- 
librium temperature of the second resistor, T2, was cal- 
culated when the temperature of the first resistor, Ti, 
was held constant. Power from only two sources was 
included: the coupling between the electrons in a re- 
sistor with the phonons at the bath temperature Tq, 



P^^ = Y,V{T^ -Tq), and photon exchange between the 



resistors and the cavity, Pr- The latter can be derived 
in the weak coupling limit where the Hamiltonian for a 
single resistor in a cavity, at position x, may be expressed 
as 



B. SIN refrigeration 

The most practical method of holding the temperature 
of the first resistor constant over a range of phonon bath 
temperatures is to remove or supply hot electrons to or 
from the normal metal by altering the bias of an SIN 
junction. Using single-electron tunnelling theory, this re- 
frigeration power can be written as"^^ 



N{E){E-eV)[fN{E-eV)~fs{E)] dE, 

(4) 

where Rt is the normal state resistance of the tunnel 
junction, fj{E) = l/[cxp( ^^ ) -|- 1] is the Fermi-Dirac 

distribution function in electrode j and N{E) is the 
superconductor density of states (DOS). Ideally there are 
no states in the superconductor energy gap. In practice 
however, it is usual to introduce a small phenomenologi- 
cal smearing parameter 7 in the DOS in order to account, 
for example, for the non-ideal superconductor or for the 
electromagnetic noise'^^''^''. The DOS is thus calculated 

as N{E) 
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straightforward to integrate Eq. (4) numerically for a 
given Tjv and Ts- 



H{x) = Q{x) ® 5Vr, 



(1) 



where SV-^os is the total voltage fluctuation across the 
resistor and Q{x) is the integral of the charge density in 
the cavity from to x. Since the coupling between the 
resistor and the cavity is linear, applying Fermi's golden 



C. SIN thermometry 

In addition, another pair of SIN junctions are needed 
on each resistor in order to measure the temperature^^. 
If a small constant current bias is maintained across the 
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junction, a measurement of the voltage allows for the 
temperature of the normal metal to be inferred. The 
SIN thermometry works on the same physical principle 
as the refrigeration and hence the current through the 
junction is'''' 
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NiE)[fNiE-eV)- fsmdE. (5) 



This allows the current bias Is to be converted to an 
equivalent voltage bias. We assume that the supcr- 
condcucting probes are thermalized with the bath at 
temperature Tq, and hence, the power introduced by 
the SIN thermometer can be calculated using Eq. (4). 
In practice this thermalization may be achieved using 
quasiparticlc traps. 



D. Photonic power transfer through the junction 
leads 



In addition to the heat conduction mechanisms dis- 
cussed above the refrigeration and thermometer probes 
also provide a supplementary channel through which un- 
wanted classical photonic heat conduction is possible. 
This excess photonic heat can occur between the resis- 
tors in the cavity or between one resistor and another hot 
element located elsewhere in the setup. We assume that 
the sample filtering and shielding is sufficiently good such 
that the latter can be disregarded and we estimate the 
former using a circuit model that consists of two resistors 
connected in series with a capacitor''^. If the tempera- 
tures of the resistors are close to each other, Ti ^ T2, the 
power through such a circuit is approximately'^'^ 
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where C is the capacitance of the series capacitor which 
is calculated assuming that each junction has a capaci- 
tance of 5 fF. This power can be reduced by fabricating 
smaller-area junctions with smaller capacitance. 



E. Quasiparticle heat transport 

The exponential decrease of quasiparticle population 
with temperature inside a superconductor implies that 
their effects should be negligible in this low-temperature 
rcgime'^^, but in principle they could make a contribution 
to the power transfer between the resistors. The heat 
power between the z*'' resistor and the quasiparticles is 
given by'^^ 
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where A is the cross-sectional area of the superconduct- 
ing line and the superconductor heat conductivity is de- 
noted by Kg. The temperature profile is governed by the 
differential equation 
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where a(T'o) « 0.98 exp (^fc^^) is the suppression fac- 
tor of the electron-phonon coupling and Eai a material 
constant for aluminium. 



F. Parameters and results 

In Fig. 2, the power into Resistor 1 from each of the 
thermal power sources in Sec. II is compared in the 
case of Tl = 100 mK and Tq = 80 mK. The temper- 
ature of Resistor 2 is scanned from 40 mK to 200 mK 
and the refrigeration bias across a single SIN junction 
at Resistor 1 is set to 14 = 0.6 A/e. The 230 resis- 
tors have a volume of 1.125 x 10^^*^ m'^ and are placed 
in a cavity of length L = 6.4 x 10"'^ m, where they 
are offset at 0.1 x L and 0.9 x L from its ends. For 
the resistors, we employ the parameters of Auo.25Pao.75 
which wc take to have a material parameter of EAu-Pd = 
3 X 10^ Wm^'^K^'^. The cavity has inductance per unit 
length of £ = 4.7 x 10^^ Hm"'^ and a capacitance per 
unit length of c = 1.3 x 10^^" Fm~^ with the fundamen- 
tal mode frequency, /i = 1 x 10^^° s^^. The normal-state 
resistance of the SIN tunnel junctions are Rt = 19 kfJ, 
the suppression of the electron-phonon coupling q;(To) = 
9.1 X 10~^, and the superconducting heat conductivity, 
Ks 1.7 X 10-1° T{x) Wm-iK-2 as used in Ref. 36. We 
use a thermometer bias current lb of 10 pA and take a 
smearing parameter 7 = 10^^. For the aluminium center 
conductor we take Eai = 0.3 x 10^ Wm-'^K"^ and an 
energy gap, A of 200 /xeV. In this setup, the normal- 
state resistance of the superconducting line with a cross 
sectional area oi A = 250 x 25 nm^ is Ri = 24 kJ7. It is 
apparent that with these parameters it is the phonon and 
photon powers that are significant, with the refrigeration 
and junction lead power being low. We may also conclude 
that the thermometer power and quasiparticle powers are 
negligible. However when the magnitude of the bias \Vb\ 
is greater than A/e, i.e., the Fermi level of the resis- 
tor is shifted above the superconductor gap, there is a 
large increase in Profr and the refrigeration power can be 
the dominant term. Since Resistor 2 has no refrigeration 
probe, this only affects Resistor 1. 



III. TEMPERATURE RESPONSE OF THE 
RESISTOR 

In the steady state, there is no net power flow at Re- 
sistors 1 and 2, and hence we have the coupled pair of 
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TABLE I. Parameters of the cavity, materials, and junctions 
which are assumed in all figures. These parameters are dis- 
cussed in more detail in Sec. IIP. Note that we always use 
identical resistors. 



Quantity 



Cavity length, L 
Inductance per unit length, I 
Capacitance per unit length, c 
Thermometer Bias Current, lb 
Junction lead series capacitance, C 
Cross-sectional area of line, A 
Smearing parameter, 7 
Resitor volume, V 
Superconductor conductivity, Ks{T) 
Superconducting gap of Al, A 
Al material parameter, Eai 
Auo.25Pao.75 parameter, Eau-pcI 
Normal state resistance of line, Ri 
Normal state NIS junction resis- 
tance, Rt 



Value 



6.4 X 10"'' m 
4.7 X 10"^ Hm"^ 
1.3 X 10""^ Fm"^ 
1 X 10~" A 
6.67 fF 
250 X 25 nm^ 
1 X 10"^ 

1.125 X 10"^° m^ 
1.7 X W'^°T Wm" 
200 fieV 

0.3 X 10^ Wm--^K- 
3 X 10'' Wm-^K-5 
24 kn. 

19 kn 
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FIG. 2. The magnitude of the power into Resistor 1 due to 
the photonic modes (upper solid line), electron-phonon cou- 
pling (dotted line), the thermometer (solid line), refrigeration 
probes (dashed line), junction leads (top dash-dot line) and 
the excited quasiparticles (bottom dash-dot line). The tem- 
perature of Resistor 2 is scanned from 40 mK to 200 mK. The 
temperature of Resistor 1 is held at 100 mK and the bath tem- 
perature is 80 mK. The resistors are offset from the ends of 
the cavity at 0.1 x L and 0.9 x L. The voltage bias across the 
SIN refrigerator is 14 — 0.6-^. All four probes are included. 
The parameters employed for the cavity are summarized in 
Table I. 
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The factors of 2 arise because there are two thermometer 
and two refrigerator probes. We do not impose a priori 
the restriction that Ti is fixed, but solve Eqs. (9) and 
(10) together to find the equilibrium values for Ti and 
T2 as a function of the refrigerator bias, Vb, on the first 
resistor. This is precisely what is done in Ref. 36 using a 



classical photonic power^-'^ for the case of two resistors in 
an impedance matched superconducting loop. The lat- 
ter setup allows the full quantum of thermal conductance 
to be achieved, i.e., the maximum possible power trans- 
fer via a single one-dimensional photonic channel. In 
Fig. 3, we observe that the temperatures resulting from 
the classical photon power used in Ref. 36 give different 
results to the quantum photonic power given by Eq. (2) 
at To 40 niK. The heating below the gap is due to the 
leakage current caused by the smearing of the density of 
states, and dissapears for 7 — > 0. 
Let us simulate the single-photon heat conduction ex- 
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FIG. 3. Resistor temperatures as functions of the bias volt- 
age on Resistor 1 (solid lines). Results calculated with the 
classical photonic power in an impedance matched loop used 
in Ref. 36 are also shown (dashed lines). With bath temper- 
ature To = 40 mK, deviations in the calculated values of T2, 
the temperature of the second resistor, are observed between 
the classical loop and quantum cavity models. The temper- 
ature is a symmetric functions of the bias. The parameters 
employed for the cavity are summarized in Table I. 



perimcnt proposed in Ref. 14 including also the refrigera- 
tion, thermometer, junction lead and quasiparticle power 
sources in addition to the phonon and photonic heat 
power. Since Eq. (10) does not depend on the refrigera- 
tion bias, with Ti held fixed, it may be solved for T2 at a 
given Ti. With this T2, Eq. (9) can then be solved to find 
the required SIN refrigeration voltage to make the equa- 
tions consistent. The results are shown in Fig. 4. The 
inclusion of the additional power sources has very little 
effect on the equilibrium value of T2. Instead the major 
consequence of incorporating the probes into the model 
is that it restricts the regime of consistent solutions for 
photonic cooling of Resistor 2. With a bath tempera- 
ture of 250 mK it is possible to find solutions only for 
Pi ^ I'^O Nevertheless, the temperature range is 

sufficient for the observation of single-photon cooling, as 
shown in Fig. 4(c). The bias voltages VJ, needed to make 
the equations consistent arc comfortably achievable in 



practice 
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FIG. 4. (a) Single-photon heating of Resistor 2 is observed 
when the temperature of Resistor 1 is swept from Ti = 40 niK 
to Ti = 400 mK for the bath temperature To = 40 mK. 
(b) Corresponding heating bias which must be applied to 
Resistor 1. (c) Single-photon cooling of Resistor 2 is ob- 
served when the temperature of Resistor 1 is swept from 
Ti = 250 mK to Ti = 170 mK for the bath temperature 
To = 250 mK. (d) Corresponding cooling bias which must 
be applied to Resistor 1. Deviations in T2 due to the effect of 
the additional thermometer, junction and quasiparticle power 
sources are indistinguishable in practice. The parameters em- 
ployed for the cavity are summarized in Table I. 



IV. PHOTON MODE DISTORTION DUE TO A 
CAPACITOR 

We first demonstrate that placing electrical compo- 
nents in a cavity changes the spatial profile of the pho- 
ton modes. This has profound effects, e.g., on the power 
transfer. We begin by simulating the system in which we 
partition our transmission line into two coupled cavities. 
We employ here the model of the transmission line as an 
infinite series of lumped LC elements. In order to divide 
this waveguide into two seperate cavities we add a capac- 
itive break of magnitude D into the center conductor of 
the coplanar waveguide as shown in Fig. 5. 
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laws at each node we obtain the eigensystem 

{-LoCouf + 2)I„ - In-i - In+i = 0, for n 7^ m,(ll) 
LoCoujf + 2 + /,„ - In^+i - I,n-i = 0. (12) 

Therefore, with no resistors in the cavities, we have an 
eigenvalue equation for the spatial dependence of the cur- 
rent, MI{xn) = LoCoutf I {xn) , whcrc the matrix M is of 
the form 



M = Mo + 



(13) 



Here Mq is the matrix without the presence of the 
capacitor which has the from of a discretized Laplacian 
and j"™ denotes the single-entry matrix with the only 
non- vanishing element equal to unity at (n,TO). We 
observe that the important factor is the ratio of the 
intrinsic capacitance Co to the coupling capacitance D. 
In the limit D — > 00 we return to the matrix for the 
unmodified cavity, and in the limit D — > 0, we have two 
independent cavities. 



This system may be solved analytically for any number 
of capacitors in arbitrary positions. In the case in which 
one capacitor is placed at L/2 we have 



A sin(fca;), 
Asm[k{L — x)]. 



for < a: < L/2 
for L/2 <x<L, 



where A is a normalisation constant and the wavevector 
k is found from 



2fc — tan — = 0. 
D \ 2 ' 



(14) 



Figure 6 shows the effect on the current profile of the 
fundamental mode when capacitors of various sizes are 
placed into a transmission line. The spatial profile is a 
sensitive function of both capacitor size and position. 
For small capacitance the current modes are vastly 
altered. As the capacitance is increased the profile tends 
back towards that of the unmodified line, as expected. 



FIG. 5. Circuit model of a transmission line cavity with a 
capacitor in the center conductor at position m. The internal 
resistance of the superconducting line may be safely neglected 
in our case. Thus the basic transmission line may be modelled 
as an infinite sequence of inductors and shunt capacitors. We 
introduce a capacitor of size L> at the position labelled by 
the index m into the central conductor. We assume that the 
inductance and capacitance per unit length, i and c, respec- 
tively, are constant along the line and hence for each inductor 
Lo = Idx and capacitor Co = cdx. 

With the capacitor placed at position index m, we 
look for solutions of the form /(x, t) = l{x) cos{u}it + (pi), 
where uJi ~ Stt/^ and 0i are the angular frequency and 
phase of the mode respectively. By applying Kirchoff's 



V. EFFECT OF THE CAPACITOR ON THE 
POWER TRANSFER 

Let us add a resistor of size R = rAx to the circuit 
of Fig. 5 at position index k. We repeat the procedure 
of applying Kirchoff's laws at this node. Due to the fi- 
nite temperature of the resistor we associate a fiuctuating 
noise voltage 5Vk{t) across it. We obtain 

- 5Vk{t) - hR - - Vk = 0, 

at 

{LoCod^ + CoRdt + 2) 4 + h-i + h+i = -CoSVkit). 

(15) 



FIG. 6. Magnitude of the current profile for tfie fundamen- 
tal mode with a coupling capacitor of size 0.1, 10, 100, and 
1000 £F (from bottom to top at x = L/2) is placed at ^. 
The corresponding values of /i are 1.96, 1.89, 1.45 and 1.07 x 
10^° Hz, respectively. The sinusoidal profile for the unmodi- 
fied transmission line is also plotted with the eigenfrequency 
fi — I X 10^° Hz. The parameters employed for the cavity 
are summarized in Table I. 



FIG. 7. (a) Normalised current profile for a coupling capac- 
itor of 0.1, 10, 100 and 1000 fF placed at ^ and resistors of 
230 SI placed at ^ and Here we assume that the resistor 
introduces no noise fluctuations into the transmission line. In 
this case, the coupling is strong enought that the values of 
/i are purely imaginary, corresponding to decaying solutions, 
(b) Normalised current profile for the same system but with 
R — 2.3 Q. The parameters employed for the cavity are sum- 
marized in Table I. 



To study the base effect on the mode profile, we first 
assume that the resistor does not introduce any noise by 
setting SV{t) = 0. In this case, Eq. (15) along with Eqs. 
(11) and (12) form the system Z(cj)/(x„) = 0. Here Z{u!) 
is an impedance matrix given by 



(w) = Mq + —J 



(16) 



For non-trivial /(x„) wc then have the condition that 
det [^(w)] = 0, yielding the eigenfrequencies of the 
cavity. This allows Z(uji) to be constructed, from which 
the current may be found. 

Figure 7(a) shows that adding a resistor of 230 VL into 
the center conductor of the transmission line changes 
the mode profile substantially. The coupling between 
the resistor and cavity is strong enough for the eigen- 
frequencies of the modes to have a dominant imaginary 
component, hence these solutions arc found to be rapidly 
decaying. In Fig. 7(b), we simulate the same system but 
with R = 2.3 0. Comparison of the two figures shows 
that the distortion of the mode is almost negligible in 
the latter case. This is likely to be the reason for the 
previously observed discrepancy between the photonic 
power transfer in the classical and quantum models in 
the strong-coupling regime^^. In order to remain in the 
weak-coupling limit, either the magnitude of the resistor 
should be small or the resistor should be positioned in 
the cavity at a point where the amplitude of the mode 
is small. In Sec. Ill and Ref. 14, the constraint that 
the resistors must remain close to the ends of the cavity 
ensures this. Otherwise, even in these simple setups it is 
possible to create vastly different spatial mode profiles 
simply by altering the position and magnitude of the 
components. 



Next we account for the voltage noise from the 
resistors. To handle 5V{t) we move into the frequency 
domain, where the current is governed by the matrix 
equation Z(lu)I{uj) = icoCoSV^uj). In our case, the 
fiuctuations are due to the Johnson-Nyquist noise, with 
the power spectrum 



Sv{uj,T) ^ 



2RhLu 



1 



exp(^)^ 



(17) 



with T being the temperature of the resistor. If a second 
resistor is placed into the line, the above procedure can 
be applied independently to both resistors in order to 
construct the impedance matrix, and thus the exchanged 
heat power may be calculated as^** 



P, 



^ J —oo 

(18) 

The results for the power transfer between two resistors 
are shown in Fig. 8 for the setup in which a capacitor in 
the center of the cavity is fiankcd by the two resistors 
at L/10 and 9L/10. The introduction of the capacitor 
reduces the power in comparison to the case of just one 
cavity. Note that with the resistors at the ends of the 
cavity, the coupling to the fundamental mode is weaker so 
that the imaginary part of the eigenfrequency no longer 
dominates for R = 230il. 



VI. TUNABLE HEAT TRANSPORT 

In this section, we show that the power transfer be- 
tween two resistors in a cavity may be tuned in situ by 
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Power transfer bet-ween two resistors, both of size 



230 f2, situated at ^5 and as a function of the magnitude 



10 10 

of the capacitor placed at ^ 



2 . The temperature of Resistor 
1 is held at 100 mK while that of Resistor 2 is at 200 mK. 
The power saturates at 2.85 fW in the limit D — >■ 00. The 
parameters employed for the cavity are summarized in Table 
I. 



placing a SQUID into the line. The transition rates for 
the absorption and emission of photons between the cav- 
ity and the resistors, given in Eq. (3) are functions of 
the effective impedance of the center conductor. This 
suggests that by inserting SQUIDs into the hne, we can 
change this effective impedance and therefore vary the 
couphng strength between the resistors in situ. 

Since we have only a few photons in the cavity we are 
comfortably in the low power regime. Following Ref. 30, 
we may then characterise a SQUID with critical current 
/co, and loop-inductance L^, as a linear inductor, pro- 
vided i^/co ^ ('J^o/27r). In this case the inductance of 
the SQUID is given by 



27r/e($) ' 



(19) 



which depends on the flux through the SQUID via the 
effective critical current, /c($) = 2/co| cos(7r$/$o)|j with 



$0 



2e 



being the flux quantum. In a real SQUID, even 



with negligible Li, there is always a slight assymmetry 
in the critical currents of the junctions Id and Ic2- As 
a consequence of this assymmetry, Ic is bounded from 
below by \Ici — Ic2\ and is therefore never identically 
zero in practice'^^. As Lj is constant for a given $ we 
may, assuming small spatial dimensions for the SQUIDs, 
transform Z{uj) — > Z^{uj) simply by replacing L with 
Lj -|- L in the circuit shown in Fig. 5 at any point in 
the line in which we place a SQUID. This is valid if the 
capacitance of the SQUID, Cj, is small enough such that 

For simplicity we take the case in which we have only 
one SQUID in the line, which also distorts the mode pro- 
file. This effect is small for most values of ^> but can be 
substantial at certain resonance points. Again we may 
use Eq. (18) to calculate the photonic power transfer. In 
Fig. 9(a), the eigenfrcquency of the cavity is shown as a 



FIG. 9. (color online) Two resistors of size 230 Q each are 
placed in the cavity at L/10 and 9I//10, between these re- 
sistors a SQUID of critical current Ico ~ 1-5 nA is placed 
at L/4. (a) Fundamental eigenfrequency of the cavity as a 
function of the flux $ with a 10 fF capacitor placed at L/2. 
(b) Power transfer between the resistors with (solid line) and 
without (dashed line) the 10 fF capacitor. The results of both 
(a) and (b) are periodic in ^ . The parameters employed for 
the cavity are summarized in Table I. 



function of the flux through the SQUID and Fig. 9(b) 
shows how the power can be modulated by the flux. 
We obtain peaks in the power exchange when on reso- 
nance. Although the maximum power may seem small, 
the power dissipated for a single photon in a very high 
quality factor cavity is of the order Pi sa . For illus- 
tration, a cavity with a quality factor of 1 x 10^ (see Ref. 
39), if we again take /i = 1 x 10^°, has Pi « 0.1 aW 
which is less than, or comparable with the powers ob- 
served in Fig. 9(b). An example of utilization would be 
to tune the cavity with the SQUID off resonance whilst 
we perform an experiment on the other half of the cav- 
ity. Afterwards, the SQUID cavity may be tuned back 
into resonance in order to reset the state of the high-Q 
cavity without the SQUID. For observation of tunable 
temperature changes in the resistor, the capacitor is not 
necessary and its removal increases the power by several 
orders of magnitude as shown in Fig. 9(b). 



VII. CONCLUSIONS 

In summary, we have used a simple circuit model in 
order to demonstrate that the photonic power between 
two components in a cavity may be tuned on and off 
using a SQUID. Furthermore, cavity experiments often 
involve components which are coupled strongly to the 
cavity modes. These component can have a substantial 
impact on the mode profiles such that they no longer bear 
any relation to the sinusoidal profile of the bare cavity. 
We have also shown that SIN thermometry does not hin- 
der the ability to observe single-photon heat conduction. 
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